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Abstract 



i -^h " In the present paper, we study an explicit effect of non-symmetry on asymptotics 

of the n-step transition probability as n — > oo for a class of non-symmetric random 
walks on the triangular lattice. Realizing the triangular lattice into M? appropriately, 
we observe that the Euclidean distance in IR 2 naturally appears in the asymptotics. 
We characterize this realization from a geometric view point of Kotani-Sunada's 
standard realization of crystal lattices. As a corollary of the main theorem, we 
prove that the transition semigroup generated by the non-symmetric random walk 
approximates the heat semigroup generated by the usual Brownian motion on R 2 . 



1 Introduction 

Let G = (V, E) be a locally finite, connected, oriented graph. Here V is the set of vertices 
and E is the set of oriented edges. For an oriented edge e G E, the origin and the terminus 



of e are denoted by o(e) and t(e), respectively. The inverse edge of e is denoted by e. A 
random walk on G is given by a non-negative valued function p on E satisfying 

= 1 for all x G V, 

e£E x 

where E x — {e G E\ o(e) = x}. Here p(e) is the probability that a particle at o(e) moves 
to t(e) along the edge e in one unit time. Then the transition probability p(n, x, y) that 
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a particle starting at x G V reaches y G V at time n is given by 



p(n,x,y) 



P(ei)K e 2)---p(e n ), 



c=(ei,e2,...,e n ) 



where the sum is taken over all paths c = (ei, e 2 , ■ ■ ■ , e n ) with t(ej) = o(ej + i), z = 1, . . . , n— 
1 and o(ei) = x, t(e n ) = y. If there exists a positive valued function my on V such that 



the random walk is said to be (my-) symmetric. 

A principal theme for random walks is to investigate the properties of p(n, x, y) as 
n — > 00. One of the most classical problems is the recurrence-transience problem which 
is related to the divergence- convergence of Y^n=iP( n i x i x )i anc ^ the l° ca l central limit 
theorem gives us a useful criterion for the divergence-convergence of the series. For this 
reason, this theme has been discussed intensively in various settings by many authors. 
See Spitzer [ID], Lawler and Limic [8] and Woess [15] for an overview of random walks. 

In particular, Kotani, Shirai and Sunada investigated long time asymptotics of p(n, x, y) 
of (my-) symmetric random walks on a crystal lattice, a covering graph of a finite graph 
whose covering transformation group is abelian. In [4J, as the precision of the local central 
limit theorem (cf. [7]), they established the asymptotic expansion 



p(n,x,y)m v (y) 1 ~ a n r/2 exp( -j- ) • (l + ai(x,y)n 1 + a 2 (x,y)n 2 H J (1.1) 



as n 00, where d(x, y) is a Euclidean distance appeared through the standard realization 
of the graph into W . In their proof, spectral theoretic arguments due to the periodicity 
of the graph and the symmetry of the random walk play crucial roles. 

Later in j2j [H] , Uchiyama and his coauthor also established the formula (II .ip for non- 
symmetric random walks on periodic graphs (i.e., crystal lattices) in the Euclidean space 
by a probabilistic approach. Their result implies that the effect of the non-symmetry on 
the coefficient a\ (x, y) highly depends not only on the underlying periodic graph but also 
on the choice of the 1-step transition probability p even if the zero mean condition (see 
condition (P2) in Section 2) is imposed. 

In view of these results, it is a meaningful problem to determine an explicit effect of 
the non-symmetry on the coefficient a\{x, y) on a specific graph. In the present paper, we 
give an answer of this problem in the case of the triangular lattice. As we will mention 
later, there exist non-symmetric random walks on the triangular lattice satisfying the 
zero mean condition (P2). On the other hand, for example, on the square lattice and the 
hexagonal lattice, we note that the zero mean condition is equivalent to the symmetry 
of the random walk since these graphs are maximal abelian coverings (see Kotani and 
Sunada [6j page 842]). In the proof of the main theorem (Theorem 12. 2p . we make use of 



p(e)m v (o(e)) = p(e)m v {t{e)), e G E, 
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the probabilistic approach as in [2J [T3] with the idea of the standard realization by Kotani 
and Sunada |U |5] . 

As a corollary of Theorem 12.21 we establish the functional analytic central limit theo- 
rem (Theorem 12 .4p that the transition semigroup generated by the non-symmetric random 
walk on the triangular lattice approximates the heat semigroup generated by the usual 
Brownian motion on M 2 . It says that the effect of the non-symmetry of the random walk 
does not appear in the appropriate space-time scaling limit. 

Throughout the present paper, O(-) stands for the Landau symbol. When the depen- 
dence of the O(-) term is significant, we specify as On(-), etc. 

2 Framework and Results 

First of all, we prepare some notations and formulate our problem. Let ei and e 2 be 
linearly independent vectors in R 2 and we set e 3 := e 2 — ei and K := ||ei ||m a + ll e 2||iR 2 - 
We define the triangular lattice G = (V, E) by 

V = {x = xiei + x 2 e 2 \ {x\, x 2 ) G Z 2 }, 

E = {(x, y) G V x V\ x - y G {±ei, ±e 2 , ±e 3 }} 

(see Figured]). We remark that G is isomorphic to the Cayley graph of the abelian group 



Figure 1: Triangular lattice 

Z 2 with a set of generators S = {±(1, 0), ±(0, 1), ±(-1, 1)} identified by V ~ Z 2 and 
e 1 ~(l,0),e 2 ~(0,l),e 3 ~(-l,l). 

Throughout the present paper, we impose the following conditions on the 1-step tran- 
sition probability p: 
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(PI): For every x G V, 

p((x, x + ex)) = a, p((x, x — e x )) = a', 

p((x,x + e 2 )) = (3', p(x,x - e 2 )) = P, 

p((x, x + e 3 )) = 7, p((x, x + e 3 )) = 7', 

where a, a', /?, 7, 7' > and a + a' + fi + ft + 7 + 7' = 1. 
(P2): Zero mean condition: 

^2p(e)e = 0. 

( P3 ) : 

r(p) := a^+^7 + 7" > 0, 
where a := a + a', (3 := (3 + f3', 7 := 7 + 7'. 

Remark 2.1 Condition (P2) zs equivalent to the following condition: 

• There exists a constant < k < | smcA i/iai a — a' = /3 — 0' — j — 7' = k. 

In tae case k = 0, condition (P2) implies thatp(e) = p(e) holds for every e G 22. Namely, 
our random walk is symmetric with my = 1 . On £/ie other hand, by a simple calculation, 
we see that there does not exist the function my satisfying ( 1 1.11) unless k = 0. Hence our 
random walk is non-symmetric in the case < k < 3. VFe can regard the constant k as 
intensity of the non-symmetry. 

We set 

M q (6):= J>(e)M> 3 , ^(^^GR^GN, (2.1) 

where (•, •) stands for the scalar product on IR 2 . Note that condition (P2) implies Mi(9) = 
0. By condition (PI), we also have 

\M q {9)\ < ( PW\\4h) PWh < K*\\0\\h, 0eR 2 ,qeN. (2.2) 

e£E 

We define the covariance matrix Q by 

(Q6,6) = M 2 {B), 9eR 2 . 

In the case of e x = ei := '(1,0) and e 2 = e 2 := '(0,1), the corresponding covariance 
matrix is easily calculated as 



Q:- 



a + 7 —7 
-7 /3 + 7 
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and hence we obtain detQ = T{p). For a general pair of two vectors ei, e 2 , we can 
decompose the covariance matrix Q as 



Q = TQ% 



(2.3) 



where T = [ei, stands for the matrix formed by column vectors ei, e 2 . It follows from 
condition (P3) and linear independence of e 1; e 2 that the covariance matrix Q is positive 
definite, i.e., (6,Q6) > A||0||jL,0 G M 2 for some A > 0. By (T2T3]) . we also observe that 



(Q 1 x, y) = (Q 1 (x 1 e 1 + x 2 e 2 ),yie 1 + y 2 e 2 } 



(2.4) 



holds for all x = x±ei + x 2 e 2 , y = yiei + y 2 e 2 G V. This means that the left-hand side of 
( 12. 4p is independent of the realization of the triangular lattice G. 
Now, we are in a position to state our main result. Let 



A(G) :-- 



^ x 1/2 
0+2 ) 
3T(p) J 



(2.5) 



1/2 



and we introduce two vectors by 

h. :- %") . , . . , . 

/3 + 7 /3 + 7 

see Figure [2]). Note that A(G) is equal to the area of the parallelogram spanned by hi 




/3+7 



/3+7 



Figure 2: hi and h 2 

and h 2 . In the case of e x = hi and e 2 = h 2 , the corresponding covariance matrix is 

Q=( 1/3 ° 
^ \ 1/3 

Realizing the triangular lattice G along these two vectors, we have 



(2.6) 
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Theorem 2.2 (1) Let us consider the case < k < 1/3. Then we have 

3 



27m-p(n,x,y) = 3A{G) exp ( - — \\y - x\\l^j 

N 

x (l + ^n-^(^)) + O^n"^ 1 ), iVGNU {0}, 

7 = 1 ^™ 



as n — )■ oo uniformly for all x, y G V , where Pj = Pj(y), j G N is a polynomial of degree 
at most 3j in the variables y\,yi and an odd or even function depending on whether j is 
odd or even. Furthermore, let us denote the leading term of J2f=i n ~^ 2 Pj {^§) > N > 2 
by ai{y — x; K)n~ l , that is, 

N 

lim n( n~^ 2 Pj ( - ) — aUy — x; n)n~ l ) =0, x, y G V. 
n^oo V ^— ' \ \/n J J 

.7=1 V 



T/ien i/ie coefficient ai(y; k) is explicitly obtained by 

ai(y; k) = af + Ko!{\y) + K 2 af\ 

where 



3T(p)i 

a? (y) = { (a? _ 2 07 + 7«)?/i + (-a? - ?7 + 2^/a)y 2 } , y = yihi + y 2 h 2 , 

fl W - 3 f i | 
1 8T( P y\ T(p) J' 

(2) Let us consider the case k = 1/3, i.e., a = (3 = 7 = 1/3, a' = f3' = 7' = 0. We set 

V k := {x = Xihi + x 2 h 2 | (xi,x 2 ) G Z 2 } + fchi, = 0, 1, 2, 

w/iere hi = 2h x - h 2 , h 2 = hi + h 2; h x ^(^2/3,0), h 2 =*(l/ v / 6, 1/^2) fsee impure Ep. 
For k,l = 0, 1, 2 ; we ta£;e x G and y £ Vi. TTien we nave 

X 1 1 TCP 1 



27m-p(n,x,y) = 9A(G) exp ( - — \\y 

N 

x (l + Y^n-^Pii^-)) +0 N (n-^), N G N U {0} 
7=1 ^ n 



as n = 3m + (/ — k) — >■ 00 uniformly for all x and y. In this case, A(G) = 1/V3 and the 
coefficient of the leading term is 

2 

a^y -x;l/3) = --. 
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Figure 3: V in the case k = 1/3 



We characterize the pair hi, h 2 through a variational problem on the crystal lattices 
based on the idea of Kotani and Sunada [HE]. 

We set ei :='(«, 0), e 2 :=*(t>i, t> 2 ). Without loss of generality, we may assume u, i>2 > 0. 
Then the energy of the (quotient graph of the) triangular lattice G is given by 

£(G) := ^p{e)\\t(e)-o{e)\\l 

e&E 

= -\^u 2 + p{vl + vl)+q{{v 1 -u) 2 + vl} . 

Minimizing £{G) with respect to (u, t>i,t>2) under the condition det((ej, e 7 -))? - =1 = A(G), 
we obtain 

U — i, f 1 = — (, f 2 = — i- 

/3 + 7 /3 + 7 

Hence we have derived ei = hi, e 2 = h 2 . 

Remark 2.3 In the case a = a' = (3 = (3' = 7 = 7' = 1 /6, £/ie random walk is said to be 
simple. In this case, we have 

, (GH _1 J = ^§, hl ^(4o), h 2 ^(-L,-L). 

VFe note the volume of the Albanese torus is equal to l/y/3 (see [U page 640], [3 Section 
10] ). The standard realization of the triangular lattice is the equilateral triangular lattice 
in M 2 each of whose edge has length a/2/3. The quantity A(G) introduced in ( 12. 5 j) can 



be regarded as a generalization of the notion of the volume of the Albanese torus to some 
non- symmetric cases. Henceforth, we call the realization of the triangular lattice G along 
hi and I12 the standard realization even if the random walk is non-symmetric. 

Next, let us consider the discrete heat equation on the triangular lattice G = (V, E) 
with initial condition Uq : V — > M: 

( (d 1 +A d )u{n,x) = 0, (n,x)eZ + xV, _ _ 

|_ u(0,x) =uq(x), x E V. 

Here diu(n, x) := u(n + 1, x) — u(n, x) and the discrete Laplacian Ad is defined by 
A d u(n, x) := ^ p(e) (u(n, o(e)) - u(n, t(e))J . 

e€E x 

The operator L :— I — is called the transition operator associated with the random 
walk on G. It is easy to see that the solution of (12.71) is given by 

u{n, x) = L n uo(x). 

We note that the solution is rewritten by 

u(n,x) = ^2p(n,x,y)u (y). 

For t > 0, let H t be the heat operator defined by 

Htf{x) = 2^/ 2 eXP (~ llX ~ 2 f R2 )f( z ) dz > f e ^(M 2 ), 

where (^^(IR 2 ) is the set of continuous functions on M? vanishing at infinity. Note that 
{H t }t>o is a semigroup whose infinitesimal generator is 

1 1 (cP_ d^_\ 

2 : ~ 2 \dxj + dx\) ' 

Now, we take the standard realization of the triangular lattice G as before. Applying 
Theorem 12.2^ we obtain the following. 

Theorem 2.4 Let t > and {5 n }^=i be a sequence of positive real numbers satisfying 
lim^oo n<5 2 = 3t. Then for every continuous function f on R 2 with compact support and 
for a sequence {x„}~ =1 in V with lim^oo 5 n x n = 16 M. 2 , we have 

lim L n (fo5 n )(x n ) = H t f(x). 

n— >oo 

This theorem is also an immediate consequence of the following result via approxima- 
tion theory due to Trotter [13] (see also Kotani [3]). 
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Theorem 2.5 (1) 



lim 

<5->0 



^A d (fo8) + l(Af)o5 



0, /eQl 2 ) 



(2) Let {<$ n }^i be a sequence as in Theorem ®. 41 

lim \\L n (f o 5 n ) - (HJ) o = 0, / G ^(M 2 ). 

n— >oo 

3 Preliminaries 

In this section, we give some basic facts for the proof of Theorem I2.2L In what follows, 
we denote h 3 := h 2 — hi, x 3 := x 2 — %i and (gfj) := (gfj) — (^7) f° r convenience. 
We begin with the explicit form of M q (9), q G N defined in (12. ip . 



Lemma 3.1 We have 
M q {6) = 



K ((e u 6)i-(e 2 ,6)i+(e 3 ,6)i) (if q is odd) , 
a(ei, 9)i + /3{e 2 , 9)" + j(e 3 , 9) q (if q is even). 



In particular, M q (9) = for every odd number q in the case where the random walk is 
symmetric. 

Noting (12. 4 p and (12. 6p . we easily see (Q~ 1 e i ,e : ,) = 3(h i( hj),z = 1,2,3. A direct 
calculation of the inner product (h*, hj), i,j = 1,2,3 yields the following lemma: 

Lemma 3.2 For any pair of linearly independent vectors ei, e 2 G R 2 ; we have 
(Q-^e,) = (Q- 1 e 2 ,e 2 ) = l±^, (Q"^, e 3 ) = 



(Q ei,e 2 ) = =7^, (Q e 2 ,e 3 ) = -— , (Q e 3) ei) 



r(p)' ^ z ' 3/ r(p)' ^ J ' x/ r(p) 

Next, we recall an elementary fact about the Fourier transform 

e -*«W> e -V=I(».*> d = 27r(detQ)- 1 / 2 e-^ Q "^^, x G R 2 . (3.1) 



Differentiating both sides of (13. ip with respect to Xi, i — 1, 2, 3, we have 

= 2ttv /z T- (detg)- 1 / 2 (- (Q^e^x^e-^ 1 ^, x G R 2 , i = 1, 2, 3. 
Repeating this argument several times, we obtain the following proposition: 
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Proposition 3.3 Let us set 



r N 

F(i u ...,i N )(x):= / (Uie^e^e-^^e-^^de, N G N, i u . . . ,i N = 1, 2, 3. 
^ / i 



T/ien we /iawe 



. . . , iiv)(x) = 2rc(V^T) N (detQ)-^e-^- 1 ^G(t 1 , i N )(x), 

where {G(ii, . . . , ik)(x) : x G M 2 ,/c = 1,...,N} is determined as the solution of the 
recursive system of the following equations starting from k — 1 to k = N: 



d 

G(h, . . .,ik)(x) = -(Q 1 e ik ,x)G(i 1 , . . .,i k -i)(x) + ^— . . .,i k -i)(x) 

G(t 1 )(x) = -(Q- 1 e tl ,x). 
Remark 3.4 G(i±, . . . ,iw)(x) is decomposed by 

f N-l 

2 



(3.2) 



G(h, . . .,i N )(x) = < 



J^G(ii,...,ijv)(2i+i)(a:) (if N is odd), 

1=0 

JV 
2 

^ G(h, . . .,in)(21)(x) (if N is even) , 



V 1=0 

where G(i\, ijy) fc = 0, 1, . . . , iV is a homogeneous polynomial of degree k in the 

variables (Q~ 1 e 1 ,x), (Q~ 1 e 2 ,x) and (Q~ 1 e 3 ,x). In particular, 

G(i 1 ,i 2 )(x) = (Q~ 1 e h ,x)(Q~ 1 e i2 ,x) - (Q~ l e h ,e i2 ), i\,ii = I, 2 , 3, 
G(i 1 ,i 2 ,i 3 )(x) = -(Q~ 1 e h ,x)(Q~ 1 e i2 ,x)(Q~ 1 e i3 ,x) 

+ (Q~ 1 e i3 ,e il )((5~ 1 e i2 ,x), ^i, i 2 , «3 = 1,2,3. 
We define the characteristic function of the 1-step transition probability p by 
<p{9) := P( e ) exp{ v /Z T(e, 9)}, 9 = (0 U 9 2 ) G R 2 . 

e£E 

We denote especially by <fi in the case of ei = ei and e 2 = e 2 . Noting that the 
characteristic function of the n-step transition probability of the random walk starting at 
the origin is equal to ip n , we have the following integral expression of p(n,x,y). 



Lemma 3.5 



p{n,x,y) 



|detT| 
(2tt) 2 

where D := [— 7r, n] 2 . 



( t T)- 1 (D) 



^(9) n exp{-^/^l(y -x,6)}d6, x,yeV, n G N, 
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Proof. Since p(n,x,y) = p(n,0,y — x), we may suppose x = 0. We also note that 
p(n, 0, y), y — y\ei + y 2 e 2 is independent of the pair of linearly independent two vectors 
ei,e2- In the case ei = e\, e2 = e 2 , the following identity is well-known: 



p(n,Q,y) 



r 2 

/ 0(0)" exp { - V=T y^y^de, y = y 1 e 1 + y 2 e 2 , n e N. (3.3) 



See Lawler and Limic [H Section 2.2.2] for details. 

For a general pair of independent two vectors e 1; e 2 , we observe 

ip(P) = <f>CT6), 9 e M 2 . (3.4) 

Then it follows from (13. 3p and ( 13. 4p that 

^)»exp{-V=T^ yi 6l 4 }de 

/ ^Tr'eye^i-v^i^y^cTr^jde 

JD i=l 



(2vr) 2 

[ ^6Texp{-V^l(y,e')}d9', y = y 1 e 1 +y 2 e 2 , n e N, 



(27T) 

|detT| 



where we performed the change of variables 9 = t T8' for the final line. This completes 
the proof. I 

Lemma 3.6 (1) In the case < k < 1/3, the characteristic function <p{9) defined on 
(*T) -1 (Z}) ^ as ^ g j \i 0V j{ n g property: = 1 holds only when 9 = 0. 

(2) In the case k = 1/3, the characteristic function <p(9) does not satisfy the above prop- 
erty. 

Proof. First, we prove (1). By (13.41) . it is sufficient to show that \<j)(9) \ = 1, 9 6 D implies 
9 = 0. We calculate the characteristic function 6(9) as 



(acos(0i) + v /Z lA;sin(^ 1 )) + (/3cos(-0 2 ) + V^k sin(-0 2 )) 
+ (7cos(^ 2 - 0i) + V^lksm(9 2 - 9J) 
=: MQ) + M0) + M0), 9 = (9i,9 2 )eD. (3.5) 

Note that the assumption < k < | implies min{a, /3, 7} > k > 0. Then we have 

<a, \<h(e)\<P, \M0)\ < 7, 9 = (9 1 ,9 2 )eD. (3.6) 
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We also observe that and |0i(#)| = a,\4>2(0)\ = j3 and 1 03 (0)1 =7 imply 0i(0) = ±a, 
2 (0) = ±f3 and 03(6*) = ±7, respectively. 

Now, we suppose |0(0)| = 1 on D. By combining (13 .5ft and ( 13. 6ft with a + /3 + 7 = 1, 
we deduce 

(0i(0), 2 (0), 03 = (a, % 7), (-«, -7), G D. 

We easily see that the first equation (0^0), 2 (0), 03(0)) = (a,(3,j) has the solution 
(6*1,6*2) = (0,0). On the other hand, there exists no solution of the second equation 
(0i(6 ) ), 02(6*), 03(6 ) )) = (— o?, — /3, — 7). Hence we conclude 6* = 0, which completes the 
proof of (1). 

The item (2) is obvious by recalling a — {3 = 7 = 1/3 and a' = (3' = 7' = 0. Actually, 
10(0)1 = 1 on D implies (9 1 ,9 2 ) = (0, 0), (f , -f ), (-f , f ). I 

Remark 3.7 In the case < k < 1/3, the random walk is aperiodic, that is, the period 
of the random walk d(p) := gcd{n G N : p(n, 0, 0) > 0} is equal to 1. It is derived from 
(1) of Lemma \3M See e.g., Spitzer [TQl P8 in Section 7] for the proof. On the other hand, 
the random walk is periodic with d(p) = 3 in the case k = 1/3. 

Before closing this section, we present an asymptotic expansion formula of the char- 
acteristic function (p which plays a crucial role in the next section. We set 

Xq (9):=(^T)- q {^Y log<p(t0), 9 G K 2 , g G N. (3.7) 



(See Bhattacharya and Ranga Rao [1, page 47].) We note that Xq{@) is a polynomial in 
the variables M 2 (9), . . . , M q (9). In particular, xi(0) = O,x 2 (0) = M 2 (0),xs(0) = M 3 (0) 
and * 4 (0) = M 4 (0) - 3M 2 (0) 2 . 

The following proposition is taken from (TJ Lemma 7.1 and Theorem 9.11]. 

Proposition 3.8 Let n G N and N G N U {0}. Then there exist positive constants 
Ci{N), C 2 (N) such that for all 9 el 2 with 



?> 1/2 < d(iv)( J2p( e )(Q^ e ) N/2 )~ N+3 

we have 



' e&E 



^(- 7 =)" - exp ( - -(Q6, 0>) ■ (b (9) + b^n' 1 ' 2 + ■■■ + b N (9)n~ N / 2 ) 

< C 2 (N) exp ( - ±(Q9, 0>) (<Q0, 9) N+3 + (Q9, 9) 3 ^ N+1 ^ ■ n^. 

Here bo(9) = 1 and bj(9), j = 1, . . . , N is written as 

b 3 {9) = {yTlf j bf\9) + (v^lf^&f ~ 2) (#) + ■ ■ • + (V^l) J+2 b { / +2 \9), (3.8) 

where bf ] (9), k = j + 2, j + 4, . . . , 3j is a polynomial in the variables Xs{@)i ■ ■ ■ > Xk{9) and 
it can be regarded as a homogeneous polynomial of degree k in the variables (e 1; 9), (e 2 , 9) 
and (e 3 , 9) . 
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Remark 3.9 In particular, 



h(9) = (v^l) 3 (^) = (v^l) 3 ( 
62(e) = (- 



24 



4 Proof of the theorems 
4.1 Proof of Theorem 1X21 

In this subsection, we prove Theorem 14. II below. Note that we easily obtain Theorem 12.21 
by combining (12. 6 j) and Lemma [3.21 with Theorem 14.11 

Theorem 4.1 (1) Lei its consider the case < n < 1/3. TTien we /iai>e 
2nn-p(n,x,y) = 3A(G) exp ( - ^(Q~\y ~ x ), V ~ x ) ) 

iV 

x (l + ^ n-^ 2 P,. (^)) + 0,v (n-^) , iV G N U {0} (4.1) 



as n —7- 00 uniformly for all x, y G V, where Pj = Pj{y), j G N is a polynomial of degree 
at most 3j in the variables 2/1,7/2 and an odd or even function depending on whether j is 
odd or even. Furthermore, let us denote the leading term of^2,- =1 n~^ 2 Pj{^^),N > 2 
by ai(y — x; k)u~ 1 . Then the coefficient ai(y; n) is explicitly obtained by 

01 (y; k) = af } + «af } (y) + K 2 aS 2) , (4.2) 

where 

aS 0) = -l + i(a(g- 1 e 1 ,ei) 2 + ^(g- 1 e 2 ,e 2 ) 2 + 7(g- 1 e3,e 3 ) 2 ), 
4\y) = \{(Q- 1 e 1 ,e 2 )(Q- 1 e 3 ,y) + (Q" 1 ^, e^Q" 1 ^, y) + (Q-WiXQ'W)}, 
4 2) = -^((g- 1 ei,e 1 ) 3 +(g- 1 e 2 ,e 2 ) 3 + (g- 1 e3,e 3 ) s 
+l((Q~ 1 ei, e 2 ) 3 + {Q- X e 2 , e 3 ) 3 - (Q^eg, 
+-((g- 1 ei, e 1 )(g- 1 ei, e 2 )(Q~ 1 e 2 , e 2 ) + (Q- 1 e 2 , e 2 )(Q~ 1 e 2 , e 3 )(g _1 e 3 , e 3 ) 
- (Q~ l e 3 ,e 3 }{Q~ 1 e 3 ,ei)(Q~ l e 1 ,e 1 ) 
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(2) Let us consider the case k = 1/3, i.e., a — {3 = 7 = 1/3, a' = /?' = 7' = 0. We set 

V k := {x = xiei + x 2 e 2 \ {x 1 ,x 2 ) e Z 2 } + kei, k = 0,1,2, (4.3) 
where ei = 2ei — e 2 , e 2 = ei + e 2 . For k,l = 0, 1, 2, we ta£;e x E Vk and y G VJ. JTien we 

2irn-p(n,x,y) = 9A(G) exp ( - — {Q~ x {y - x), y - x) ) 

x ( 1 + Z) n " i/2 ^(^7^)) +°iv(n-^), iVG NU{0}(4.4) 
i=i ^ n 

as n = 3m + (I — k) — > 00 uniformly for all x and y. In this case, the coefficient of the 
leading term is 

2 

ai(y-x; 1/3) = --. 

Proof. As we mentioned in the proof of Lemma I3.5[ we may suppose x = throughout 
the proof. 

First of all, we prove (1). By recalling that the covariance matrix Q is positive definite 
and Lemma I3.6[ we can choose a positive constant R sufficiently small such that both 

r, := sup{|^(6»)| : 9 G (T^D), ||d|| R3 > R} < 1, (4.5) 

and 

\<p(0)\ <exp(-i(Q0,0», P\W*<R (4-6) 

hold. See e.g., pllj P7 in Section 7] and also Shiga [9] Lemma 6.15] for the proof of ( 14.61) . 
By Lemma 13.51 we have 

2tt 



'(*T)-i(C) 

^(detQ) 1 ^ / J-^-Ye^^de, (4.7) 

J y/n{*T)- x (D) 



2tt 



where we performed the change of variables \/nd' = 9 and used (12. 3p and (12. 5 \ for the 
second line. We take a positive constant r sufficiently small such that 

2N+3 N 

r < mm{C 1 (N)K-^+T\^+^,R}, (4.8) 

and divide the range of the above integration ^/n( t T)~ l (D) into three parts according as 
2 < rn 1 / 6 ; rn 1 / 6 < ||0||r2 < Ry/n; \\9\\r2 > Ry/n. Then we can write as 



/ K 

J^T)-^(D) V 



9 V " -V=T<2/,-?=>J/) 

e v ^'V^'dO 



I{n){y) + Ji(n)(y) + J 2 (n)(y) + J 3 (n)(y) + J 4 (n)(2/), 
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where 

J(n)(y) 



N N 



J!(n)(2/) := / {J^Y -e-^^Yb^n-i/Me-^'^dO, 

„ AT 

J 2 (n)(y) := - / e-^^fVfo^^n-^e-^'^^, 

Hn){y) := / ^(-^)V^'^ ) ^, 

J 4 (n)(y) := / ipl-^Ye-^^dO. 



Our first task is to estimate the error terms Ji(n)(y), J2(n)(y), J^{n){y) and J±(n)(y). 
By using (14. 5 p and (14. 6p . we have 

sup \Mn){y)\ < [ f] n d9 = (27r) 2 |detT| • mf, 

y&v J^T)-^(D) 

and 

sup | J 3 (n)(y) | < / e--^ Q9) d6< [ e-* m l*d8 

Thus both J 3 (n)(y) and J±{n){y) converge to zero as n — > oo exponentially fast uniformly 
for all y G V. Similarly, we have 

N 

sup\J 2 (n)(y)\<e-^ nl/4 Yn- j / 2 ( / e"^ 11 * 2 \bj(6)\de) , 

and since each has polynomial growth, we also see that J2{n)(y) converges to zero 
as n — > oo exponentially fast uniformly for all y G V. 
By (T22D, we have 

p(e)(Q~ 1 e, e) N l 2 < K N X~ N / 2 , NeNU {0}. 

e£E 

Furthermore recalling (14.81) . we observe that \\0\\m? < rn 1 ^ 6 implies 
(Q6,6) 1/2 < K\\e\y 

N N 1 
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Thus we may apply Proposition 13. 8[ and we obtain 

Bup|Ji(n)(y)| < C 2 {N)n-^ [ e -|W> ((Q9,9) N+3 + (Q9,9) 3 ^)d9 

„ , N+l . 

< N (n- — ). 
Now, we calculate the principal terms 

I^iy) := n' 3 ' 2 [ ^^^^^"^i 1 ^, y e V, j = 0, 1, . . . , N. 



'M 2 

It follows directly from (13. ip that 



I {n)(y) = 27r(detg)- 1 / 2 exp ( - l-(Q^y,y)). 



Applying Proposition 13. 3[ we obtain 

■M 3 (9) 



Ji(n)(y) = (v^l) 3 ™- 1 / 2 / (^^)e~^ e ^e~^ I{y >^ ) d9 

Jr2 V 6 / 

= 27m- 1 / 2 (detQ)- 1 /2 e xp ( - ^-(Q^y^Px^) 



2n \/n J 
where 

Pi(y) ■= /t(V g Il)6 1, - G(2, 2, 2)(y) + G(3, 3, 3)(y)). 
Besides, it follows from Remark 13.41 that 

G(z, 2, z) (^=) = -n" 3/2 (Q- 1 ^, y} 3 + Zn- 1 ' 2 {Q~ 1 e l , e t ) (Q~ l e h y) , z = 1, 2, 3. 

Thus the explicit form of Pi J is given by 

Pi{^=) = ^- 3/2 {(Q- 1 ei,y) 3 -(Q- 1 e 2 , 2 /) 3 +(Q- 1 e 3 ,y) 3 } 

+ ^- 1 / 2 {(g- 1 e 1 ,e 2 )(g- 1 e 3 ,y) + (Q- 1 e 2 ,e 3 )(g- 1 e 1 ,y) 

+ (Q- 1 e 3 ,e 1 )(g- 1 e 2 ,y)}. (4.9) 

It follows from Remark 13.91 that 

I 2 (n)(y) = n 2 ' 2 [ (v^l) 4 ( ^ - *) e 4W> e -^4>| 
y R2 V 24 8 / 

./m>2 V 72 / 
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Then by repeating the same argument as above, we have 

I 2 (n)(y) = 27m- 2 / 2 (detQ)- 1 / 2 exp ( - ^(Q^V, v)) (p£\n)(^=) + P«( n )( * )), 
where 

P?\n)(y) = ^=^{aG(M, l,l)(y) + ?G(2, 2,2,2)(y) 
+7^(3, 3, 3, 3) (y)} 

^{a 2 G(l, 1, 1, l)(y) + £ 2 G(2, 2, 2, 2)(y) 



+7 2 G(3, 3, 3, 3)(y) + 2af3G(l, 1, 2, 2)(y) 
+2^G(2, 2, 3, 3)(y) + 2^aG(l, 1, 3, 3)(y)}, 

A (2) (n)(y) = ^ ( ^T 1)12 { E G(i, i, i, i, i, i) (y) - 2G(1, 1, 1, 2, 2, 2) (y) 



72 

i=l 



+2G(1, 1,1,3, 3, 3)(y) - 2G(2, 2, 2, 3, 3, 3)(y)}. 
Besides, it follows from (13.21) that 



G(i,i,i,j,j,j)(o)(y) = -9(Q 1 e i ,e i )(Q 1 e i ,e j )(Q 1 e j ,e j )-6(Q ^e,-) 3 , 
G(i,i,j,j)( )(y) = (Q~ 1 ei,e i )<Q~ 1 e i ,e i ) + 2<Q- 1 ei,e i ) a , i, j = 1,2,3. 



By combining these identities with Lemma [3.21 it holds 

3(S 2 (Q" 1 e 1 , ei > 2 + PiQ- 1 **, e 2 ) 2 + f^"^, e 3 ) 2 ) 
+2S / S((g~ 1 e 1 , ei ) (Q- 1 e 2 , e 2 ) + 2<g- 1 e 1 , e 2 ) 2 ) 
+27a((g- 1 e 1 , e x ) (Q^eg, e 3 > + 2<g- 1 e 1 , e 3 > 2 ) 
+2^7((g- 1 e 2 , e 2 ) (Q^eg, e 3 ) + 2(Q~ 1 e 2 , e 3 ) 2 ) = 8. (4.10) 

Thus the constant term of P 2 and P 2 are obtained by 

- 1 + 1 (a<Q- x ei, ei ) 2 + ^Q- 1 ^, e 2 > 2 + ^(Q-^, e 3 ) 2 ) (4.11) 

and 



((g^ei, ei ) 3 + (Q~ l e 2 , e 2 ) 3 + (g^eg, e 3 ) ; 
+~ ((g-^x, e 2 ) 3 + (g-^a, e 3 > 3 - (g^eg, ei ) 3 ) 

+ ^ ((g^ei, e 1 )(g- 1 ei, e 2 )(g- 1 e 2 , e 2 ) + (Q -1 e 2 , e 2 )(g^ 1 e 2 , e 3 )(g- 1 e 3 , e 3 ) 

-(g- 1 e 3 ,e 3 )(g- 1 e 3 ,e 1 )(g- 1 e 1 ,e 1 )), (4.12) 
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respectively. Summarizing (14 . 9 j) . (I4.10p . (14. lip and (I4.12p . we conclude that the explicit 
form of the coefficient of the leading term is given by (14.21) . 

For general j — 1, 2, . . . , N, by virtue of (13.81) . we may write Ij(n)(y) as 

h(n){y) = j^lf\n){y) 

k=l 

■= ^{^ /2 (v^l) J+2fe / $ + *\e)e-?WVe-^ v >3* ) d0}. 
k=i " /r2 

Applying Proposition 13.31 again, we obtain 



lf\n){y) = 2n(^iy +2k n-^(detQ)-^e W ( - {Q~ l y, y))Pj j+2k) (-^=) , (4.13) 
where 

J2 (^iy +2k - 2l p^M (if j is odd), 



/ 3-1 



J 1 fc 



if +2fc) fo) := 



z=o 
i+fe 



J2(^l) j+2k ~ 2l P^ 2 lM (if j is even) , 
v 1=0 

and each Pjj^k-2i(y) ls a homogeneous polynomial of degree (j + 2k — 21) in the variables 
y u y 2 . Noting (j + 2A; - 2/) < 3j and (v^l)^ 2 ^^!)^ 2 ^ 2 ' G R for any j, k, I, we see 
that 

^(y):=E(v^l) J+2fc ^ +2fc) (y), J = l,2,...,iV (4.14) 

k=l 

is a real valued polynomial of degree at most 3j in the variables y±, y 2 and it is an odd or 
even function depending on whether j is odd or even. 
Then by fT4TT3l) and ( Oil) , we obtain 

I 3 {n){y) = 2 7 m^/ 2 (detQ)- 1 / 2 exp ( - ^(Q' 1 y,y))P j (^y=). (4.15) 



Here we mention that the term n~ J / 2 Pj(^) on the right-hand side of (I4.15P is regarded 
as a polynomial of degree at most 2j in the variable n -1 . 

Plugging the above all arguments into (14.71) . we finally obtain the desired asymptotic 
expansion formula (14. ip . This completes the proof of (1). 

Next, we prove (2). For simplicity, we only give the proof in the case k = I = 0. 
(In other cases, the proof goes through in a very similar way with a slight modification.) 
Let ]?(3) := p * p * p and we denote by and Q(3) the characteristic function and the 
covariance matrix associated with the 1-step probability distribution p( 3 ), respectively. 
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We define \\ (0), geNin the same way as in in (13.TB with ip replaced by cpp). Then we 
easily see Qr$\ = 3Q and Xq (y) — $Xq(p)i Q e N. 

Let G = (V,E) be an enlarged triangular lattice defined by V := V and E : = 
{(x,y) G V x Vo| x — y G {±e 1; ±e 2 , ±(e 2 — ex)}} (see Figure EJ). Note that 

xiei + x 2 e 2 = yiei + y 2 e 2 = ziei + ^ 2 e 2 

implies 

s)-Ks)- (s) =s (*) := (-i !)(:)■ (4 - i6) 

Now, we consider a random walk on G whose 1-step transition probability distribution 
is given by prg\ . We denote by (m, x, the m-step transition probability of the random 
walk. Noting P(3)(m, x, y) := p(3m,x,y), x,y G V and combining Lemma 1331 with (14.161) . 
we have 

27r(3m)-p(3m,0,i/) = ^H\det(TS)\ f m {9) m e~^ 1 ^d9 

= M^)( det Q)l/2| det5 | 

x [ m (^J) m e -^y^)dB, (4.17) 



Here, we observe the random walk is aperiodic because of p( 3 \(l, 0, 0) = 2/9 > 0. By 
virtue of Spitzer [TUl P8 in Section 7], we see |<^(3)(#)| = 1,9 G {'(TS')} -1 ^) implies 
= 0. Thus we may follow the proof of (1), and by 1 14. 17ft . we obtain 



27r(3m) .p(3m,0,y) = ^^(detg) 1/2 |detS'| V /,-(m)(y) + O^m"^ 1 ) (4.18) 

27T Z — ' 



3=0 

as m — > oo, where 



Ij(m)(y) := mT j/2 I b j (9)e~^ Q ^ e ' e) e~ V ^ I{y '^ ) d9, y E V,j = 0,1, . . . , N 



■ — ' Sn\ 

and bj(9) is defined in the same way as (13. 8 j) with x g (6 ) ) replaced by Xg (#)• 

Furthermore, repeating the calculation of Ij(n)(y) and recalling Q(3) = 3Q, we obtain 

5(m)(y) = 27r(3 2 detg)- 1 /2 exp { - ^-.(Q-'y, y)}^-^). (4.19) 

We finally obtain the desired asymptotic expansion formula (14. 4p by combining (I4.19P 
with (I4.18P and noting | det | = 3. This completes the proof. I 
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4.2 Proof of Theorem I2T41 

First, let us suppose that < k < 1/3. Recall that the sequence {^n}^^ satisfies nS 2 — > 3t 
as n — > oo. Since the area of the fundamental domain in M 2 with respect to Z 2 -action of 
5 n G is S 2 n A{G), 



H t f(x) = lim J-Ve-All^n-MI| 2/(W5 2 A(G) 

n->oo z — ' 
ydV 



lim 



E 



e -|rll x ™-fll, ; 



' /(*»!/)• 



Theorem 12.21 implies that, for any e > 0, there exists iV G N such that for all n > N, 
\n8t — 3t\ < t and 



(3A{G) _i 



27m 



2 JI*« y\\ R2 - p (n,x n ,y) f(8 n y) 



< 



Since / has a compact support, by the definition of A(G) in (I2.5P and the property of S n , 
3A(G)£ J2\f(S n y)\ =^J2\f^y)\€A(G) 



27m 



< 



n yeV 

3e 



: 2tt (3t- |n<J2-3i|) 



£|/(<S^)|<fr4(G) 



3e 
~4vrt 



£|/(kl/)|<&4(G)- / |/(*)|cfe] + 



3e 
4vrt 



\f(z)\dz 



for some C > 0. Then we conclude Theorem 12.41 in the case < k < 1/3. 

Next, let us consider the case k — 1/3. For n G N, let m n G N be the quotient of n 
divided by 3 and k n G {0, 1, 2} be its reminder. Namely, m n and k n satisfy n = 3m n + k n . 
For the sequence {s„}™ =1 C V satisfying 5 n x n — > x G M 2 as n — > oo, let {i n }nLi be a 
sequence of 0, 1, 2 satisfying x n G Vi n , a subset of given in (14.31) . By the periodicity of 
the random walk, we have 



y)f{S n x n ) 



y^ x ^y)-^^- e ^ l]Xn ~ yll] 



f(S n y) 



+ 



2im 



E 



/((J„y)3^(G), 
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where V3 = Vo, V4 = V\. By using (2) of Theorem 12.21 and the previous argument for 
< k < 1/3, the first term of the right-hand side converges to as n — > 0. Since the area 
of the fundamental domain in M 2 with respect to Z 2 -action of 5 n Vk n +i n is 35^A(G), 

lim -A- e-^-<>f(5 n y)38 2 n A(G) = -L / exp ( - ^ " f ll » a )/( g )cfe. 

This completes the proof of Theorem 12.41 I 



4.3 Proof of Theorem 1231 

First, let us prove (1), i.e., the convergence of the infinitesimal generator. By using 
Taylor's theorem with respect to the coordinate (xi,X2) — x\hi + X2h2, we have 

^A d (fo8)(x) = ±£ Ke) Wx) - f(5t(e))) 

e£E x 

e<EE x i=l,2 4 

i,j=l,2 2 ^ 

= -^E^( fe )E^ e ) e 

8=1,2 * eEE x 

-\ E E^)^)-^^)-^^:^) 

i,j=l,2e£E x 1 j 

+0{5). (4.20) 

Thanks to the zero mean condition (P2), the first term of the right-hand side vanishes. 
Since the covariance matrix Q with respect to the standard realization is given by (12. 6ft . 
the second term of the right-hand side is — |(A/)(&c). Hence we conclude (1). 

The item (2) is deduced from (1) by applying Trotter's approximation theory [13] (see 
also [3]). This completes the proof. I 



5 Examples 

Example 5.1 (Kotani-Sunada J3j, Example 2) We consider the simple random walk. 
Namely, we consider the case a = a' = (3 = (3' = 7 = 7' = 1/6. As we mentioned in 
Remark 12.31 the standard realization of the triangular lattice is the equilateral triangular 
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lattice in M 2 each of whose edge has length a/2/3. Then we have 
27m-p(n,x,y) = \/3exp ( - ^-\\y - x\\l^ 

x(l-k 1 + '" + ^(^)) +0 "("^)' N ^ 3 
as n — y oo uniformly for all x, y G V. 

Example 5.2 ( |12j ) We consider the case 

a = p = 7 = (1/6) +e, a' = p' = 7' = (1/6) - e (0 < e < 1/6). 

In this case, the standard realization of the triangular lattice is same as the above example. 
However, the corresponding random walk is non-symmetric with k = 2s. Then we have 

2im-p{n,x,y) = v^exp ( - —\\y - x\\l^ ■ {l + ( - - - Qe^jn' 1 

+ '" + ^(^f)} + ^"^)- JVS3 

as n — > 00 uniformly for all x = xihi + x 2 h 2 , y = yihi + y 2 h 2 E V. 
Example 5.3 We consider the case 

a = 1/4, a' = 1/12, p = 1/3, p' = 1/6, 7 = 1/6, 7=0 (k = 1/6). 
In this case, 

2 , 2V2 . , /2V2 \ . , ( 1 1 
7TT' 7TT' hl= (7iT' j' h2= iTl'Tf 

(see Figure HJ) . Then we have 

. , 6 / 3 ll2 

2irn-p(n,x,y) = -j= exp ^ - — - x\y 

■ f 719 6 

1+ - 



x 



1331 121 
1 „ (y — x s 



fl»(J^)] +0w ( n -*P), iV>3 



n 7V/2- 

as n — ?• 00 uniformly for all 2 = xihi + X2h2, ?/ = yih.\ + y2fr-2 G V. 

Acknowledgement. The authors would like to thank Professors Atsushi Katsuda, Mo- 
toko Kotani, Kazumasa Kuwada and Yukio Nagahata for helpful discussions and com- 
ments. 
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Figure 4: Example 5.3 
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